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We discuss the optimal detection strategy for a stochastic background of gravita- 
tional waves in the case n detectors are available. In literature so far, only two cases 
have been considered: 2- and n-point correlators. We generalize these analysises 
to m-point correlators (with m < n) built out of the n detector signals, obtaining 
the result that the optimal choice is to combine 2-point correlators. Correlating n 
detectors in this optimal way will improve the (suitably defined) signal-to-noise ratio 
with respect to the n = 2 case by a factor equal to the fourth root of n(n — l)/2. 
Finally we give an estimation of how this could improve the sensitivity for a network 
of multi-mode spherical antennas. 
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I. CORRELATION OF TWO DETECTORS 

As it is well known the sensitivity to a stochastic background signal can be greatly 
enhanced by correlating the output of two detectors. To show how this works it is useful to 
consider the cross correlation Syi [2] between two detector outputs S\ and 5*2, defined by 

/T/2 i-T/2 poo 

dt / dt' S 1 (t)S 2 {t')Q{t -t')= dfS*(f)S 2 (f)Q(f) , (1) 
-T/2 J -T/2 J-oo 

where the filter function Q(t) has been introduced. The cross correlation Si 2 depends only 
on the time difference t — t' as stationarity in both the signal and the noise is assumed. 
In the last equality the Fourier transform of the signal and the limit T — > oo have been 
taken. For any finite T, Su is made of the sum of statistically independent random variable 
involving Si(f) and S^/'), which are correlated only over a frequency range \f — f'\ < 1/T. 
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Thus, as Si 2 is the product of random variables, it is a random variable itself and it can be 
approximated by a Gaussian variable by virtue of the central limit theorem, even in the case 
of narrow band detectors, provided that T is much larger than the inverse of the bandwidth. 
The same will be true in the case the product of more than two random variables, that will 
be considered later. 

The outputs of two detectors can be split as Si )2 = si >2 + iV lj2 , being Sj the physical 
signal and iVj the noise. The signal-to-noise ratio for the correlation of the 2 detectors at 
our disposal (this redundant notation will be useful later, where m-point correlators out of 
n detectors will be considered) is given by 

[SNR(2|2)f = <M = W = fog* , (2) 

where (S i2 ) and o"i 2 are respectively the average and the square root of the variance of the 
cross correlation. We have adopted the convention which makes the signal-to-noise ratio 
proportional to the metric perturbation h, so that in our notation SNR oc h, as in j^j], 
differently from |2j where SNR oc h 2 . To obtain the last equality in (j2J) we have made 
the basic assumptions that we will never drop throughout this paper: both the signal and 
the noise are Gaussian, they are statistically independent, stationary and with zero mean, 
Ni ^> Si and finally the noises of different detectors are completely uncorrelated. 

The filter function Q(t) appearing in (JJJ can be freely chosen in order to maximise the 
signal-to-noise ratio. The best choice is obtained in the standard way by imposing the 
functional variation of (0) with respect to Q(t) equal to zero and solving it for Q(t). To 
write down the explicit form of the filter function it is necessary to introduce some further 
quantity. The signal can be usefully written as 

/oo p 
dfi / dnjtAifi, n^^-^F^) , (3) 
-oo J 

where Ha is the Fourier transform of the metric perturbation with polarisation A and F A is 
the pattern function of the detectors, which encodes the information on its angular sensitiv- 
ity, Xi the position of the z-th detector and the wave arrival direction. Given the stochastic 
nature of the signal, the 2-point correlator (ensemble average of the Fourier components) of 
the metric perturbation can be parameterised as 

Wu ni)/i(/ 2 , na)> = 8{f x + f 2 )^5 2 (fli, fla) \s h {h) , (4) 
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where the spectral function Sh has been introduced. Analogously a noise spectral function 
Sjv,i for the z-th detector can be defined through 

(N^N^h)) = s^ih + f2)\s N Ah) • (5) 

The filter function which maximises the signal-to-noise ratio is 

where the overlap function F has been introduced. Its definition involves the relative distance 
and orientation of the two detectors 

r(f t ,x ab ) = l.f d 2 nJ2F^\n)F^\n)e 2m ^-^ , (7) 

being F^"^ (Q) the pattern function of the detector at site a, b for a wave coming from 
direction Cl. Inserting the optimal filter function © in and in (J2J) the explicit form of 
the signal-to-noise ratio for the correlation of two detectors is obtained 

s 2 h (f) ^ 1/4 



/ poo 

SNR(2,2)= T / dfT 2 (f,x 12 )— - 



(8) 



which gains in the case of two identical detectors with respect to the single detector case, 
as it is well known, a factor roughly equal to (TA/) 1//4 multiplied by the overlap function, 
being T the experiment time and Af the bandwidth. 

II. CORRELATION OF N DETECTORS 

One now might ask what can be gained by the correlation of several such detectors. A 
partial answer is obtained by generalizing (JSJ) to the case of 2n detectors (the number of 
detectors must be even for the correlator not to vanish) Q] 



SNR(2n|2n) (^•■■ 5 *» 



((Nl..N%)f 2 

T i/±( f d * f ^ (lir=i Sh(fi)T(fj, x i)n+i )f + perm 

\J J U2=l^N,i(fi)SN,i+n(—fi) 



l/4n (9) 



where in our notation SNR(z|j) is the signal-to-noise ratio given by i-point corre 
out of j detectors. To obtain Q the explicit form of the optimal filter function 



ators taken 
, which is 
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determined up to an arbitrary constant, 

Q(f2, ■ ■ ■ , hn) OC 

(10) 

S h (-f n+1 )T(-f n+1 , a?i, w +i) [IIL2 s h(fi)T(fj, x i:n+l )5(f n+i + fj)] + perm 

5jV,l(— fn+l) Yl i= 2 ^N,i{fi) 

has been used. The permutations come from all the different pairings of 2n signals: as 
the detector output is Gaussian its n-point correlator can be computed from the product of 
2-point ones. In SNR(2n|2n) we indicated only the terms with leading behaviour in T for 
large T, which are (2n — 1)!!. Eq. (JHJ) can be rewritten schematically as 

n 

[SNR(2n|2n)] 4n = JJ [SNR(i, n + i\2)f + perm , (11) 

8=1 

i.e. it is the sum of permutations of products of 2-point correlators. 

We now show that there exists a better way to treat data obtained from 2n detectors, as 
out of 2n detectors, 2m-correlators can be considered, for any m < n. For m = 1 we can 
follow the analysis of Q or Q and consider all the possible pairs taken out of 2n detectors. 
For each detector pair a mean value and a variance can be defined as usual 

= (Sij) = S 2 afj = (Si) - Sfj , (12) 

where the optimal filter function has been normalized so to make the theoretical mean 
(S^) = S 2 equal for every pair. A SNR(z, j\2) of the type (jSJ) can thus be assigned to each 
pair 

[SNR(z,j|2)] 2 = ^. (13) 

Oij 

The best way to gather the information from all the pairings is to take a weighted average 
with weights \j 

S 2 = %± X f J , (14) 



whose variance is 



*l = (si) - (s 2 y = E ^ A ^ 2 , 



(J2i<j x v) 
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which is justified by large noise approximation we are using, that allows to neglect non 
diagonal terms like cr^-er^ (for {i,j} ^ {A;, I}) compared to afy The signal-to-noise ratio 
obtained by combining in pairs the 2n detector outputs in this way is given by 

[SNR(2 | 2 „ )f= ^>! = (^yj . (15) 



a 



S 2 ^i<j ij ij 

The best signal-to-noise ratio is obtained by choosing Ay oc cr^ 2 (which correspond to weigh- 
ing less the more noisy data) and it is 

[SNR(2|n)] 4 = |f = E [ SNR (^il 2 )] 4 > ( 16 ) 

i<j l 3 i<j 

where we have dropped the unnecessary hypotheses of the number of detectors being even. 
The optimal signal-to-noise ratio is thus given by the sum of terms like (jHJ) (to the fourth 
power); note that we recover the time dependence of (J^J: SNR(2|n) oc T 1 / 4 For n 
detectors with equal noise level, data collection time and overlap functions, we have 

SNR(2|n) oc [n{n - 1)] 1/4 . (17) 

We now generalize the analysis of the combination of 2-points correlators to the case of 
2m-point correlators. Analogously to (jl4j) we can define a linear combination S2 m of the 
(n)!/[(2m)!(n- 2m)!] 2m-point correlators S^...^ that is possible to build out of n detectors. 
Defining a signal-to-noise ratio of the type © 

[SNR(zi . . .i2m\2m)\ = = , 

as a natural generalization of ()13|h we are led to consider the combination of the 2m- 
correlators analogous to (fTlj) 

V A <? V A 2 n 2 

q _ Z^n<...<i 2m /V tl-»2m' J 2m 2 _ Z^i 1 <...<i 2 m il-i2m U il-i2m 

Z^h<...<i 2m /X U...l2m l2^» 1 <„.<i Sm A *l-i2mJ 

(with z'fc e {1 . . . 2n}) so that the signal-to-noise ratio for 2m-correlators can be written, for 
the optimal choice of weights \i 1 ...i 2m ^ a i x ...i 2 ■> as 

[SNR(2m|2n)] 4m = = £ (^O! = ^ [SNRfa . . . * 2m |2m)] 4m . (18) 

tl<...<t 2m ^l-^m h<...<i 2m 

Each of the terms in the sum in the most rhs of (fTBj) is on its own the sum of (2m — 1)!! 
terms as shown in 
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For equal noises, observation times and overlap functions the scaling of the signal-to-noise 
ratio with respect to the number of detectors n, and with the order of the correlator 2m, is 
given by 

i 

(19) 

where the first factor comes from the number of contribution in each SNR(ii . . . i 2 m,|2m) and 
the binomial coefficient from the possible choices of 2m-ple out of n detectors. 

For any fixed n the maximum is obtained always for m = 1 implying that the optimal 
signal-to-noise ratio is obtained by combining the detectors in pairs as in (j!6|) . In particular 
for a network made of a large number of detectors the signal-to-noise is expected to scale 
with the square root of the number of detectors as in (|T7jl . 



SNR(2m|n) oc 



(2m- 1)!! x 



n 



2m 



III. STOCHASTIC BACKGROUND 



To parameterise conveniently the detector sensitivity to a stochastic background it is 
useful to introduce the normalized spectral energy density of gravitational waves Q gw (f) 
defined as follows 

<W/) = -%^, (20) 
pc d In/ 

where p c = {?>h\ m Hl) / (87rG n) is the critical energy density of the Universe (h W0 H Q is the 
Hubble constant and Hq = 100Mpc _1 km/s) and p gw is the Fourier transform of the energy 
density in gravitational wave. In terms of the spectral function Sh introduced in (HJ) p gw can 
be written as 

~d^T = 2G~J ShU) ■ (21) 
Using this formula it is possible to rewrite (jH)) in terms of Q gw 



[SNR(2|2)f = V2T 



th'LM \ /"°^r 2 (/,x 1 , 2 )^( / ) 



'-loo-'-'o 
4tt 2 



L/0 f e £>N,lSN,2 



1/2 

(22) 



which can be used in fjT7)j) to express the SNR(2|n) as a function of Q gw . The SNR necessary 
to claim detection can be computed once a false alarm rate a and a false dismissal rate 
1 — 7 are specified (7 is called detection rate), according to the Neymann-Pearson detection 



criterion (see j^j for the definition of false alarm and false dismissal rates). A natural choice 
is a = 5% and 7 = 95%. 

Even if it is useful in principle to correlate as many detectors as possible, in practice 
the number of high sensitivity detectors available with non negligible overlap function is 
not too large. For instance in j^l it is shown that with the current light interferometers 
it is possible to detect a stochastic background of gravitational waves provided that their 
normalized energy density satisfies the bound 

^f°' 5% >6.5-10- 6 , (23) 

for 7 = 95% and a = 5%, constant Q gw and an observation period of 4 months. This bound 
is obtained by correlating the two LIGO's, VIRGO and GEO600, but a numerically similar 
one is obtained by correlating just the two LIGO's. 

From the phenomenological point of view it is interesting to note that the total gravita- 
tional wave energy stored in a stochastic background cannot exceed the bound 

tt 9W < 6 • 10" 6 , (24) 

surprisingly numerically close to ()23j) . otherwise the Universe would expand too rapidly 
in the epoch of primordial nucleosynthesis, thus spoiling the beautiful agreement between 
theory and observation of the primordial abundance of light elements. The nucleosynthesis 
epoch corresponds to ten seconds after the Big Bang or to a temperature of the order of few 
MeV. 

This limit does not apply to a background of gravitational waves produced after the 
nucleosynthesis epoch. There exist indeed astrophysical sources which may produce a con- 
tinuous stochastic signal in the phenomenologically interesting frequency range. For instance 
rapidly rotating young neutron stars could be the source of gravitational radiation with an 
amplitude of h\ m VL gw ~ 10~ 8 for the frequency range 0.7 — 1 kHz |4j. 



A system which is able of producing gravitational 



sented by cusps and kinks in cosmic string network 5]. For some range of the parameters 
entering the description of the cosmic strings the energy of gravitational wave Q gw could be 
as high as 10~ 6 and it could be produced either before or after the nucleosynthesis epoch. 



wave to an observable level is repre- 
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IV. NETWORK OF ANTENNAS 

We have shown in sec. |H] that with n detectors available, the best strategy to detect 
a stochastic background is to correlate pairs of them, and in case noise levels and overlap 
functions have equal values the signal-to-noise ratio increases as n 1//2 for large n. 

Let us now turn our attention to a case in which it is important to correlate a large 
number of detectors to increase the detector sensitivity. 

An interesting case can be realized by multi-mode detectors like spherical antennas [f| as 
they have five quadrupolar modes which couple to a gravitational wave with generic incident 
direction. The correlation of such modes among several antennas can be considered, thus 
increasing the number of effective detectors available. 

Anyway it has to be considered that modes of the same sphere cannot be correlated 
among themselves, as their noises are correlated and most of the mode pairs have negligible 
values of the overlap functions. Fig. shows the overlap reduction function for different 
pairs of modes. Denoting by m the integer number labelling different quadrupolar modes 
(—2 < m < 2), the five overlap functions in the figure are obtained by correlating the m = 
mode of a sphere and the five modes of a second sphere located at 100 Km (the quantisation 
axes have been chosen parallel to each other in order to maximise the sum of the overlap 
functions). 

This makes the statistics a bit different than in (|17|). so that for a set of n spherical 
antennas the analogous of (fTHjl is 

5 

[SNR sp/iere (2|n)] 4 = [SNR(*,m;j,m'|2)] 4 , (25) 

i<j<n m,m'=l 

where the indices i,j run over different detectors and m, m' run over the quadrupolar modes. 
If all the modes are equally noisy it is clear from fig.^that every mode of the first sphere can 
be effectively correlated with just one or two modes of the second sphere, as it has almost 
vanishing overlap with the others. The situation does not improve by choosing a different 
orientations for the quantization axis of different spheres, as an example for an angle of 45 
degrees shows in fig. El 

At the moment the only operating spherical antenna is miniGRAIL (and another one will 
soon start working, the Brazilian Gravitational Wave Detector Mario Schenberg j^), thus 
to obtain some numbers let us consider the present noise spectral function of miniGRAIL 
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Figure 1: Overlap reduction function, as a function of the frequency, between the m = mode of 
a sphere at lat. 52.2° N, long. 4.5° E, (quantisation axis chosen along the local meridian) and the 



five quadrupolar modes (labelled 20, lc, Is, 2c, 2s after 



) of a second sphere located at 100 Km 



(lat. 52.2° N, long. 6.0° E, quantization axis along the local meridian). 



.71 



whose diameter is 68 cm. It can presently reach a strain sensitivity h c of about h c = 
y/Sh ~ lO^^Hz^ 1 / 2 , it has a resonant frequency of 2.9kHz and a bandwidth of about 230 
Hz. For this figures one obtains, for a single pair of detectors 



/ 2 95%,5% 
' t 100 si ffu> 



1 1 4ir z 



df 



r 2 (/,x 12 ; 



-1/2 
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(26) 



for an observation time of 1 year. Using a set of identical spheres eq. (|25|) can be explicited 

as 



[SNR sp/iere (2|n)] 4 = 2T 



(27) 



4vr 2 7 Jo / 6 ^(/) 
where Sjv(/) is the common noise spectral function and r(/, Xjj) is the overlap function for 
the detector pair z— j, which is understood to depend also on m and w! . The situation can be 
considerably improved by using larger spheres, with a consequent lower resonant frequency. 
We can estimate for instance that slightly improving the sensitivity to h c ~ 10~ 21 ifz~ 1//2 , a 
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resonant frequency at, say, 300 Hz (and bandwidth 100 Hz) one can reach 

"1/2 

(28) 

which is obtained by inverting (|27j) for a constant Q gw , thus allowing to obtain Q gw ~ 3 • 10 -5 
in the experiment bandwidth for a set of n = 10 detectors. We note that it is important to 
have a not too high resonant frequency for detector correlation, as overlap functions go to 
zero at / > 1/nL, being L the detector separation. This still far from light interferometers, 
see eq.()23|) and the phenomenological bound given by eq.(|24|). but the effect of the multiple 
correlation is quantitatively important, so it is not excluded that once higher sensitivity will 
be achieved the correlation effect will be crucial for detection. 

The mechanism of sensitivity enhancement actually will not work if a sphere is correlated 
with an interferometer as only one mode of the sphere can effectively be correlated with an 
interferometer. This can be seen in fig. El which shows the overlap function between VIRGO 
and the five quadrupolar modes of a hypothetical sphere placed in Rome. A sphere like one 
with the characteristics leading to 1)28)1. has a narrower bandwidth than an interferome- 
ter, but similar sensitivity in its bandwidth, so correlating a sphere with an interferometer 
would lead to a result equal to (f2~Kj) but for the factor in square brackets, as one cannot 
take advantage of the correlation of several modes of the same sphere: correlation with an 
interferometer would just add one more single-mode detector. 

V. CONCLUSION 

We have analyzed the utility of considering multiple detector correlators to detect a 
stochastic background of gravitational waves. The main result of the paper is the demon- 
stration that the best way to correlate the outputs of different detectors is in pairs, no 
matter how large is the number of detectors, instead of taking m-correlators with m > 2. 
Finally as a potentially interesting application of this result we applied this strategy to a 
set of identical spheres, showing that correlation of several pairs of detectors is important 
in increasing the sensitivity. 
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Figure 2: Overlap reduction function, as a function of the frequency, between the m = mode of 
a sphere at lat. 52.2° N, long. 4.5° E, (quantization axis along the local meridian) and the five 
quadrupolar modes (labelled 20, lc, Is, 2c, 2s after p) of a second sphere located at 100 Km 
(lat. 52.2° N, long. 6.0° E, quantization axis at 45° with respect to the local meridian). 
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